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A  New  Design  of  Differential  Space-Time  Block  Code 
Allowing  Symbol- Wise  Decoding 

Yu  Chang,  Yingbo  Hua,  Fellow,  IEEE,  and  Brian  M.  Sadler,  Eellow,  IEEE 


Abstract — For  four  (or  more)  transmitters,  a  new  design  of 
differential  space-time  block  code  allowing  symbol-wise  decoding 
is  presented  in  this  letter.  The  new  design  not  only  has  the 
minimum  (symbol-wise)  decoding  complexity  as  that  by  Yuen, 
Guan  and  Tjhung  (YGT)  but  also  yields  a  lower  error  rate.  While 
the  YGT  code  uses  a  specially  designed  symbol  constellation,  the 
new  code  uses  a  conventional  QAM  with  a  rotation.  At  a  high 
data  rate  snch  as  3bps/Hz,  the  new  design  with  symbol-wise 
decoding  complexity  can  even  yield  a  lower  error  rate  than  the 
code  by  Zhu  and  Jafarkhani  that  has  the  pair-wise  decoding 
complexity. 

Index  Terms —  Coding-gain,  differential  space-time  code, 
symbol-wise  decoding. 


I.  Introduction 

IFFERENTIAL  space-time  block  codes  (DSTBC)  are 
useful  for  wireless  communications  with  multiple  trans¬ 
mitting  antennas.  With  DSTBC,  the  channel  state  information 
(CSI)  is  not  required  either  at  the  transmitters  or  at  the 
receivers,  which  is  important  for  applications  where  the  CSI 
changes  too  fast  to  be  estimated  and  utilized.  The  design  of 
DSTBC  has  attracted  the  attention  of  many  researchers  in 
recent  years.  For  two  transmitters,  the  design  of  DSTBC  is 
well  established  because  of  the  existence  of  full  rate  complex 
orthogonal  code.  But  for  more  than  two  transmitters,  the 
design  of  DSTBC  is  still  an  active  area  of  research.  For 
practical  use,  there  is  a  strong  interest  to  reduce  the  decoding 
complexity  of  DSTBC  with  as  little  loss  of  coding  gain  as 
possible.  Consider  a  system  of  four  transmitters  with  a  full 
diversity  DSTBC  operating  at  the  data  rate  (i.e.,  spectral 
efficiency)  r  bps/Hz.  A  conventional  decoder  requires  a  size- 
2'*’'  search,  a  pair-wise  decoder  requires  a  size-2^’'  search,  and 
a  symbol-wise  decoder  needs  only  a  size-2’'  search.  Although 
a  pair-wise  decoder  reduces  the  complexity  significantly  from 
a  conventional  decoder,  the  complexity  difference  between  a 
pair-wise  decoder  and  a  symbol-wise  decoder  is  still  large.  For 
four  transmitters,  while  many  designs  of  pair-wise  decodable 
DSTBC  are  now  available  in  the  literature  such  as  the  codes  by 
Zhu  and  Jafarkhani  (ZJ)  [1]  and  Calderbank  et  al  [2],  the  code 
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by  Yuen,  Guan  and  Tjhung  (YGT)  [3]  seems  the  only  reported 
DSTBC  that  is  of  full  rate,  full  diversity,  allows  symbol-wise 
decoding,  and  suffers  no  major  loss  of  coding  gain  from  the 
pair-wise  decodable  counterparts. 

In  this  letter,  we  present  a  new  design  of  DSTBC  that 
allows  symbol-wise  decoding.  The  new  design  differs  from 
the  YGT  code  in  the  following  ways.  The  YGT  code  forces 
the  symbol  constellation  to  be  constructed  in  such  a  way  that 
a  symbol-wise  decoder  becomes  near  optimal  with  respect  to 
the  given  encoder.  The  resulting  symbol  constellation  consists 
of  symbols  falling  on  the  real  and  imaginary  axes,  and  the 
constellation  spread  along  one  axis  is  wider  than  that  along 
the  other.  See  Fig.  3  in  [3].  The  new  design  utilizes  a  conven¬ 
tional  QAM  constellation  and  applies  a  rotation  to  optimize 
a  “pseudo  coding  gain  (PCD)"  governed  by  the  encoder.  The 
decoder  of  the  new  design  is  forced  to  be  symbol-wise  in 
an  approximate  maximum  likelihood  fashion.  Although  the 
decoder  of  the  new  design  is  not  optimal  with  respect  to 
its  given  encoder,  the  higher  PCD  of  the  encoder  allows 
the  loss  of  the  suboptimal  decoder  to  be  more  than  enough 
compensated.  As  shown  by  simulation,  the  error  rate  of  the 
new  design  is  significantly  smaller  than  that  of  the  YGT  code 
at  moderately  high  SNR.  While  the  idea  of  the  new  design 
is  simple  which  innovatively  combines  those  of  YGT  and  ZJ, 
the  new  design  also  shows  the  highest  coding  gain  among 
all  known  DSTBC  that  allow  symbol-wise  decoding  for  four 
transmitters. 

Next,  we  show  the  details  of  the  new  design  as  well  as  its 
connections  to  the  YGT  and  ZJ  codes.  The  simulation  results 
are  illustrated  in  Section  III. 


II.  The  New  Design 

WE  start  with  a  constellation  set  S'  =  {si,i  = 
1,2,...,M}  in  the  complex  plane  where  M  =  2’’  . 
We  assume  that  there  are  21  transmitting  antennas  where  I  is 
an  integer-power  of  2. 

If  the  actual  number  a  of  transmitting  antennas  satisfies 
29  <  a  <  29+1  ,  one  should  treat  the  number  of  transmitting 
antennas  to  be  29+i  .  The  difference  29+i  —  a  can  be  treated 
as  the  number  of  channels  that  have  zero  gain  among  total 
29+1  nominal  channels.  The  transmit  diversity  factor  remains 
to  be  a. 

During  the  time  block  of  21  time  slots,  we  collect  2m 
symbols  randomly  from  S  .  This  process  should  be  governed 
by  an  outer  layer  code  that  is  not  discussed  here.  These 
symbols  are  denoted  by  =  {sk,i,i  =  1,2,  ...,2m}  where 
m  <  I  . 
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A.  First  Step  of  Encoding 

We  construct  a  2m  x  1  symbol  vector  as  follows: 


Xfc  = 

[  ^k,l  ^k,2  f 

Xfc.l  = 

[  ^k,l  ^k,m  ] 

Xfc,2  = 

[  ^k,2m  ] 

^k,i  — 

and  i  =  1,2,  ...,m  .  The  mapping  from  Sk  to  x^  is  the  same 
as  a  part  of  the  YGT  encoder.  See  the  entries  in  the  quasi- 
orthogonal  matrix  in  Equation  (8)  in  [3]. 

A  key  property  of  the  YGT  code  is  also  governed  by  a 
quasi-orthogonal  code  first  shown  in  [4].  It  is  shown  in  [5] 
that  the  quasi-orthogonal  code  shown  in  [4]  is  a  member  of  a 
large  family  that  all  hinge  on  the  2x2  orthogonal  code  first 
shown  in  [6]. 


Therefore,  we  can  also  write 

ak,i  = 

After  a  substitution  of  the  encoder  described  earlier,  it  follows 
that 

CTfc.l  = 

Oik, 2  = 

It  is  important  to  note  that  for  i  =  1,2  is  a 

linear  superposition  of  such  terms:  (3?(sfe_t)  ±  for 

t  =  1,2,..., 2m  .  This  property  is  critical  for  symbol-wise 
decoding  as  shown  later. 


B.  Second  Step  of  Encoding 

We  now  construct  a  21x21  code  matrix  that  is  partitioned 
into  four  I  x  I  submatrices: 

(Xfc),,=Yfc.i;(Xfc),2  =  Yfc,2; 

(Afc)2  2  =  (Afc)2  2  = 

where  Y^^  =  (Afc,i  -f  Afc,2)/2  and  Yfc,2  =  (Afc,i  -  Afc,2)/2. 
The  matrices  A^,  i  and  Akg  are  encoded  differentially  as 
follows: 

Afc  1  = - Mfc  iAfc_i  1 

y/ak-is 

Afc  2  =  —  Mt-  2Afc_i  2 

^Oik-1,2 

where  ak-i,i  =  det{Al^_i  ,^Ak-i^i)  ,  is  an  orthogonal 

matrix  M(zf,  i)  of  the  symbol  vector  z^  i  =  x^  i  -f  x^  2,  and 
Mfc,2  is  an  orthogonal  matrix  M(zfc  2)  of  the  symbol  vector 
Zfc,2  =  Xfc  i  —  Xfc  2-  At  the  beginning  of  each  frame  that  is 
comprised  of  multiple  blocks,  we  set  Ao,i  =  Ao,2  =  I  with  I 
being  the  identity  matrix.  This  part  of  the  encoder  resembles 
the  differential  encoder  by  ZJ:  see  (28)-(31)  in  [1]. 


C.  A  Combined  Effect 

We  limit  M(z)  to  be  a  linear  matrix  function  of  the  vector  z 
and  of  the  property:  M^(z)M(z)  =  ||z|pl  with  ||z|p  being  the 
squared  norm  of  z.  For  I  =  2,  M(z)  can  be  a  full-rate  Alamouti 

matrix  M(z)  =  ,  for  which  I  =  m  =  2  .In  this 

—  Z2  Zi  J 

case,  there  are  four  transmitters  as  21  =  4,  and  the  rate  of 
symbols  per  time  slot  is  one  (i.e.,  full  rate).  When  I  >  2  , 
there  are  more  than  four  antennas  as  21  >  4,  and  we  have  to 
use  fractional-rate  (i.e.,  I  >  m)  orthogonal  matrices  as  shown 
in  [7]  and  the  references  therein. 

It  is  easy  to  verify  that 


A^iAfc,!  =  =  ||zfc,ifl 

A^2Afc.2  =  M^2Mfc.2  =  llz^fl 


D.  The  Channel  Model 

It  is  our  convention  that  each  column  of  the  code  matrix 
Xfc  contains  symbols  transmitted  from  one  antenna,  and  each 
row  of  the  code  matrix  X^  contains  symbols  transmitted  during 
one  time  slot.  We  assume  a  single  receiver  while  the  multiple- 
receivers  case  can  be  easily  handled  as  mentioned  later.  Then, 
the  received  signal  vector  during  the  time  block  can  be 
written  as:  =  X^h  -|-  n„  where  h  is  the  channel  response 

vector  from  21  transmitters  to  a  single  receiver,  and  is  the 
noise  vector  that  is  assumed  to  be  white,  Gaussian,  and  have 
the  variance  .  We  assume  that  h  is  quasi-static,  i.e.,  the 
change  of  h  from  one  time  block  to  the  next  is  negligible. 
However,  h  does  not  need  to  be  constant  during  a  whole  time 
frame  in  order  for  the  differential  code  to  be  used,  although 
our  simulation  shown  later  assumes  a  constant  h  during  each 
frame  for  the  purpose  of  comparison  to  other  codes.  We  also 
assume  that  all  elements  in  h  are  independent  and  identically 
distributed  complex  Gaussian  random  variables  that  change 
independently  from  frame  to  frame. 

E.  Pseudo  Coding  Gain 

We  partition  the  received  vector  r^,  into  two  I  x  1  subvectors: 
Tfc  1  and  Tfc  2  ■  Similarly,  h  is  partitioned  into:  hi  and  h2  ,and 
Hfc  into  nfc  i  and  Then,  using  the  definition  of  the  code 
matrix  X^,  it  is  easy  to  verify  that 

^k,i  =  ffc,! +  rfc.2  =  Afc,i(hi -f-h2) -f  nfc,i -f  nfc,2 
^k,2  =  ffc.i  -  rfc,2  =  Afc,2(hi  -  h2) -f  -  nfc,2 

Note  that  q^,  1  and  2  ^re  sufficient  statistics  of  the  un¬ 
known  symbols  as  they  are  a  unitary  (up  to  the  factor  \/2 
)  transformation  of  the  original  observations  1  and  Yk^2  ■  If 
the  channel  response  vector  h  is  available  at  the  receiver,  the 
matrix  Ak  =  diag{Ak^i,Akg)  serves  as  an  equivalent  code 
matrix.  Hence,  the  minimum  value  of  the  following: 

(det  (^AAfAfc^^^ 

=  (^det  (^AAj^^^Ak^i'j  det  (^AA^^^kg'jY^ 

=  {Aak,iAak,2)^  (3) 
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with  respect  to  the  differential  symbol  set  AS'^  =  {sk^i  — 
Sk,j  7^  0;  *  =  1)  2, 2to;  j  =  1,2,...,2to}  ,  would  be  the 
coding  gain  of  an  optimal  coherent  decoder.  Here,  the  operator 
A  implies  the  use  of  differential  symbols.  For  example,  given 
that  Ak  is  a  function  of  the  symbol  set  Sk  ,  AA^  is  the  (same) 
function  of  the  differential  symbol  set  ASk  ■  For  this  reason, 
we  define  a  pseudo  coding  gain  as 


G 

=  minfcmin  ^  {Aak,iAak,2)^ 

1  <  j  <  m 
1  <  i  <  m 


=  minAseAS 


[5i(As)]^ 


[9(As)]^ 


for  which  the  symbol  constellation  is  normalized  to  have 
unit  variance.  The  pseudo  coding  gain  becomes  the  exact 
coding  gain  of  a  differential  decoder  under  some  conditions 
as  discussed  later. 


F.  The  Decoder 


To  derive  a  symbol-wise  differential  decoder,  we  need  to 
exploit  the  recursion  of  A^.  i  and  Afc_2  in  q^.  i  and  2  >  which 
leads  to 

Qfc.i  =  Mfc.iQfc-i,!  +  Wfc,i 

qfc,2  =  ,  ^  Mfc,2qfc-i,2 +Wfc,2  (4) 


where 


Wfc,i  = -  Mfc4(nfc_i  1  +  nfc_i^2)  +  (nfc,i  +  T^k,2) 

Wfc,2  = -  Mi- 2(nt,_i  1  —  nfc_i^2)  +  (nfc,i  —  nk^2) 

^C(k-l,2 

It  can  be  shown  that  the  composite  noise  vectors  1  and 
Wfc_2  are  each  white,  Gaussian  and  have  variances 


^2 

=  20-2  1 

^2 

^fc,2 

=  20-2  1 

Ctk-1,1  ) 
Oik, 2  \ 
ctk-1,2  ) 


(5) 


respectively.  Furthermore,  Wk,i  and  Wfc^2  are  independent  of 
each  other  under  the  Gaussian  assumption  of  . 

If  ctfc.i  =  Oik, 2  =  const  hence  ^  =  o'fe  2  =  const,  the 
optimal  (i.e.,  maximum  likelihood)  differential  decoder  of  the 
block  is  given  by 


ming^  J  =  minsJJi  +  J2) 


where 

~  Mfc^aq^_2  j|| 

^Ok-l,i 

for  7  =  1, 2  .  In  this  case,  the  pseudo  coding  gain  G  becomes 
the  exact  coding  gain  [8].  If  ak,i  =  ak,2  ^  const  (i.e.,  they 
are  dependent  on  the  symbols  in  the  block),  the  above 
detection  is  near  optimal,  and  G  may  be  treated  as  a  near 
coding  gain.  In  fact,  under  the  YGT  constellation  [3],  ak,i  = 


Fig.  1.  A  4-QAM  constellation  with  a  rotation  angle. 


Ok,2  7^  const  and  the  above  detection  becomes  equivalent  to 
the  YGT  decoder. 

If  Ok,i  ^  Ok,2  ,  then  min5j,(Ji  +  J2)  is  only  an  approxi¬ 
mation  of  the  maximum  likelihood  detection,  and  G  becomes 
only  a  guide  of  the  actual  coding  gain.  Unfortunately,  the 
maximum  of  G  over  all  possible  constellations  is  achieved 
when  ak,i  7^  ak,2  ■  This  implies  the  existence  of  tradeoff 
between  the  value  of  G  and  the  optimality  of  min5^(Ji  +  J2) 

.  The  best  tradeoff  should  yield  the  highest  actual  coding  gain. 
The  YGT  code  in  [3]  chooses  one  extreme  of  the  tradeoff 
where  min5^(Ji  +  J2)  is  forced  to  be  optimal  conditionally 
upon  a  chosen  set  of  constellations.  This  penalizes  the  value 
of  G.  The  best  value  of  G  under  ak,i  =  ak,2  7^  const  can  be 
shown  to  be  1.33  for  2  bps/Hz  and  0.49  for  3  bps/Hz,  which 
is  consistent  with  Fig.  2  in  [3].  But  for  the  4-QAM  and  8- 
QAM  constellations  (see  Figs.  1  and  2)  with  the  rotation  angle 
6  =  0.5ton“^(0.5)  «  13.28°  ,  the  G  values  can  be  shown  to 
be  1.79  for  2  bps/Hz  and  0.60  for  3  bps/Hz.  These  values  are 
significantly  larger  than  those  for  the  YGT  constellation.  The 
rotation  angle  13.28°  is  optimal  in  maximizing  G  under  the 
QAM  constellation,  which  is  also  shown  in  [9].  For  example, 
referring  to  Fig.  2,  one  can  verify  that  at  0  =  0.5fan“^(0.5)  , 


G 


(3?(Asab))^  —  (9(Asab))^| 

(Sft(AsBD))^  —  (Q(Asbd))^ 

(3?(Ascb))"  -  (9(Ascd))" 
0.60 


The  square-shaped  8-QAM  constellation  shown  in  Fig.  2  is 
not  the  only  8-QAM  constellation  that  achieves  G  =  0.60  .  In 
fact,  a  rectangular-shaped  8-QAM  that  has  8  symbols  on  the 
corners  of  3  consecutive  squares  also  has  the  same  G  value 
under  the  rotation  angle  13.28°  .  But  the  corresponding  ak,i 
and  ak,2  are  statistically  less  close  to  each  other,  which  makes 
minsj.  ( Ji  -f  J2)  less  optimal.  Consequently,  the  actual  coding 
gain  using  the  rectangular-shaped  8-QAM  is  not  as  good  as 
using  the  square-shaped  8-QAM,  which  has  been  confirmed 
by  our  simulation. 

Note  that  in  general,  a/c,i  and  ak,2  are  random  variables 
under  different  k.  One  way  to  characterize  the  difference 
between  them  is  to  use  statistical  distributions.  It  can  be 
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Fig.  2.  An  8-QAM  constellation  with  a  rotation  angle.  Here,  A,  B,  C  and 
D  indicate  the  locations  of  four  symbols.  The  difference  between  the  two 
symbols  at  A  and  B  is  denoted  by  As^s- 


shown  that  under  a  random  selection  of  symbols  of  zero 
mean,  i  and  2  have  the  same  distribution.  But  what 
is  more  important  for  the  optimality  of  min5^(Ji  +  J2) 
is  the  difference  between  cr^  ^  and  ^  ■  This  difference 
would  be  small  if  ak.i  and  ak,2  do  not  change  much  as  k 
changes.  From  the  expressions  of  1  and  ak,2  in  terms  of 
the  original  symbols,  we  can  see  that  if  m  is  large  then  a^.i 
and  ak,2  should  be  stable  and  hence  cr^  ^  and  cr^  2  should 
be  close  to  each  other.  In  our  simulation,  we  have  observed 
that  even  for  m  =  2  ,  the  distributions  of  ^  and  cr^  2 
tend  to  concentrate  toward  a  common  constant.  This  makes 
min5j,(Ji  +  J2)  not  a  bad  approximation  of  the  maximum 
likelihood  detection.  However,  it  is  worth  stressing  that  the 
optimality  min5^(Ji  +  J2)  always  depends  on  the  symbol 
constellation. 

It  is  critical  to  show  that  regardless  of  its  optimality, 
min5j,(Ji  +  J2)  can  be  always  decomposed  into  independent 
symbol-wise  decoders.  Recall 

Ji  =  \\qk,^  -  ^ 

y/C(k-l,i 


for  i  =  1,2.  Then,  it  follows  that  for  i  =  1,2: 


J^ 


Hk.iV- 

Q!fe-l,z 


2 


2 


(6) 


From  (1)  and  (2),  we  know  that  al/l  is  a  linear  superposition 
of  such  terms  of  ±  Q{sk,t))^  for  f  =  1, 2, ...,  2m.  We 

also  know  that  Mj,  ^  is  also  a  linear  superposition  of  terms  of 
individual  symbols  Sk^t  for  f  =  1, 2, ...,  2m  .  Therefore,  J1+J2 
can  be  straightforwardly  decomposed  into  a  sum  of  symbol- 
wise  cost  functions,  and  each  cost  function  can  be  minimized 
independently. 

If  there  are  multiple  receiving  antennas,  there  is  a  cost 
function  Ji  +  J2  for  each  receiving  antenna,  and  the  sum 
of  all  these  cost  functions  is  then  equivalent  to  (exactly  if 
afc.i  =  Oik, 2  =  const  ,  nearly  if  ak,i  =  ctk,2  ^  const  , 
or  approximately  if  ak,i  Q!fe,2  )  the  likelihood  function  of 


Fig.  3.  Block  error  rates  for  2  bps/Hz.  In  descending  order  of  error  rates  at 
SNR=25dB,  shown  here  are  the  YGT  code,  clairvoyant  YGT,  proposed  code, 
its  clairvoyant  version,  the  ZJ  code,  clairvoyant  ZJ,  and  the  SP(2)  code.  The 
difference  between  the  proposed  code  and  the  ZJ  code  is  almost  negligible. 
Note  that  SP(2)  uses  a  search  of  size-256,  the  ZJ  code  uses  a  search  of  size- 16, 
and  the  proposed  code  and  the  YGT  code  use  a  search  of  size-4. 


the  total  received  signals.  The  total  cost  can  be  similarly  de¬ 
composed  into  symbol-wise  cost  functions.  The  total  diversity 
becomes  the  product  of  the  number  of  transmitters  and  the 
number  of  receivers. 


III.  Simulation 

We  assume  four  transmitting  antennas  and  one  receiving 
antenna.  Each  frame  consists  of  64  blocks  and  each  block 
has  4  symbols.  The  channel  response  (fading)  factors  are 
constant  during  each  frame  but  vary  from  frame  to  frame  as 
independent  complex  Gaussian  random  variables.  The  noise  is 
white  Gaussian.  The  SNR  of  the  received  signals  is  computed 
hy  E  {\\Xkh\y)  /E  {\\nkr). 

Fig.  3  shows  the  block  error  rates  of  four  codes  versus  SNR 
for  2  bps/Hz.  The  first  is  the  YGT  code  and  its  clairvoyant 
version.  The  clairvoyant  version  is  such  that  the  symbols  at 
block  fc  —  1  are  perfectly  known  when  the  symbols  at  block 
k  are  detected.  The  second  is  the  proposed  code  and  its 
clairvoyant  version.  Both  the  YGT  code  and  the  proposed  code 
are  decoded  in  the  symbol-wise  fashion.  The  constellation 
used  for  the  proposed  code  is  the  rotated  4-QAM  as  shown 
in  Fig.  1  with  9  =  13.28°  .  The  third  code  is  the  ZJ  code 
and  its  clairvoyant  version.  This  code  uses  pair-wise  decoder. 
The  fourth  code  is  the  SP(2)  code  that  is  decoded  by  a 
program  provided  by  authors  of  [10].  For  the  2bps/Hz  case, 
the  SP(2)  code  requires  a  size-256  search,  the  ZJ  code  needs 
a  size-16  search,  and  both  the  YGT  code  and  the  proposed 
code  need  only  a  size-4  search.  But  the  difference  among  the 
coding  gains  of  the  four  codes,  even  with  different  decoding 
complexities,  is  not  very  large.  It  is  important  to  note  however 
that  the  proposed  code  has  a  consistently  better  coding  gain 
than  the  YGT  code. 

We  now  consider  the  quaternion  code  by  Calderbank  et  al 
shown  in  [2].  This  code  requires  size-16  search  for  2bps/Hz. 
The  block  error  rate  of  the  code  is  not  available  in  [2],  but 
is  shown  in  Fig.  4  for  2bps/Hz  -  courtesy  of  S.  Das  of  [2]. 
We  see  that  when  SNR=20dB,  the  block  error  rate  is  about 
9  X  10“"^  which  is  better  than  all  above  codes  except  the  SP(2). 
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Fig.  4.  Bit  error  rate  and  block  error  rate  of  the  quaternion  code  by 
Calderbank  et  al  for  2  bps/Hz  -  Courtesy  of  S.  Das  of  [2]. 


Fig.  5.  Block  error  rates  for  3  bps/Hz.  The  proposed  code  and  the  YGT 
code  use  a  search  of  size-8,  and  the  ZJ  code  uses  a  search  of  size-64.  In  this 
case,  the  proposed  code  outperforms  both  the  YGT  code  and  the  ZJ  code  by 
a  significant  margin  at  high  SNR. 

But  when  SNR=25dB,  the  block  error  rate  is  about  4  x  10“® 
which  is  nearly  the  same  as  the  YGT  code  but  worse  than  all 
other  codes.  It  appears  that  the  slope  of  the  error  curve  of  the 
quaternion  code  is  not  as  steep  as  others. 

Fig.  5  compares  the  block  error  rates  of  the  proposed  code, 
the  YGT  code  and  the  ZJ  code  for  the  3bps/Hz  case.  For 
3bps/Hz,  both  the  proposed  code  and  the  YGT  code  use  a 
size-8  search,  and  the  ZJ  code  requires  a  size-64  search.  The 
constellation  used  for  the  proposed  code  is  shown  in  Fig.  2 
with  the  rotation  angle  9  =  13.28°  .  The  constellation  used  for 
the  ZJ  code  is  the  same  as  Fig.  2  but  with  6  =  tt/4  -  courtesy 
of  Y.  Zhu  of  [1].  The  ZJ  code  with  9  =  13.28°  did  not  perform 
as  well  as  with  0  =  7r/4  .  The  8-PSK-constellation  was  also 
considered  for  the  ZJ  code.  But  it  yielded  poorer  results.  We 


see  that  the  block  error  rate  of  the  proposed  code  is  about 
one  third  of  that  of  the  YGT  code  at  SNR=25dB.  We  also  see 
that  the  ZJ  code  becomes  even  worse  than  the  YGT  code  at 
SNR=25dB.  This  is  somehow  different  from  the  result  shown 
in  Fig.  5  in  [3]. 

Finally,  we  note  that  only  at  1  bps/Hz,  the  ZJ  decoder  was 
shown  to  be  the  exact  maximum  likelihood  decoder.  For  this 
reason,  it  seems  no  surprise  that  at  a  higher  data  rate  such  as 
3  bps/Hz,  the  ZJ  code  can  perform  worse  than  the  proposed 
code  and  the  YGT  code.  A  bit  error  rate  comparison  of  the  ZJ 
code  and  the  proposed  code  was  also  conducted  at  3  bps/Hz, 
which  yielded  a  similar  conclusion. 

IV.  Conclusion 

Symbol-wise  decodable  DSTBC  (differential  space-time 
block  codes)  are  important  for  practical  applications  in  wire¬ 
less  communications  with  multiple  transmitters.  For  four  trans¬ 
mitters,  the  YGT  code  [3]  is  symbol-wise  decodable  and  has 
a  coding  gain  close  to  the  ZJ  code  [1]  that  requires  pair-wise 
decoding  [1].  The  code  proposed  in  this  paper  is  also  symbol- 
wise  decodable  and  outperforms  the  YGT  code  by  a  significant 
margin.  The  new  code  results  from  a  fusion  of  the  YGT  code 
and  the  ZJ  code. 
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